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Abstract 

We establish an algorithm that produces a new solution to the Ein- 
stein field equations, with an anisotropic matter distribution, from a 
given seed isotropic solution. The new solution is expressed in terms of 
integrals of known functions, and the integration can be completed in 
principle. The applicability of this technique is demonstrated by gen- 
erating anisotropic isothermal spheres and anisotropic constant den- 
sity Schwarzschild spheres. Both of these solutions are expressed in 
closed form in terms of elementary functions, and this facilitates phys- 
ical analysis. 
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1 Introduction 



Solutions to the Einstein field equations with anisotropic matter, where the 
radial component of the pressure is not the same as the tangential component, 
have been studied by a variety of investigators in recent years. Such solu- 
tions are important in physics as they may be used to model anisotropic stars 
where a static spherically symmetric interior is matched to the Schwarzschild 
exterior model. These include the treatments of Dev and Gleiser 0| , Her- 
rera et al. 0, J|, Ivanov j^, Mak and Harko ^ 2| and Sharma and Mukherjee 
0. As pointed out in these analyses anisotropy cannot be neglected when 
studying the critical mass and redshift of highly compact bodies, and may 
be an intrinsic component of boson stars and strange stars. 

In generating exact models describing anisotropic matter distributions ad 
hoc assumptions are normally made about the forms of the matter variables 
or gravitational potentials. In an attempt to systematically generate exact 
solutions Maharaj and Maartens Q proposed an algorithm in which the field 
equations are written as a first order system of differential equations. The 
energy density and radial pressure are then chosen on physical grounds; the 
remaining variables then follow from the field equations. The Maharaj and 
Maartens algori thm has been used by Gokhroo and Mehra jl^, Chaisi and 



Maharaj |lll.ll2| amongst others to produce physically reasonable anisotropic 
stars. However these models have the feature that the anisotropy factor 
is nonzero in the interior of the star and an isotropic solution cannot be 
regained. The objective of this paper is to provide a new algorithm that may 
be used to systematically produce anisotropic solutions to the Einstein field 
equations from a given isotropic solution. This solution generating algorithm 
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has the desirable property of yielding a new exact anisotropic model which 
includes an isotropic solution, unlike the earlier algorithm of Maharaj and 
Maartens. The new solutions obtained are members of a one parameter 
family containing the isotropic case. We illustrate the utility of the algorithm 
with two examples that are physically relevant, the isothermal sphere and 
the constant density Schwarzschild sphere. 

The principal objective of this paper is to show that we can generate 
anisotropic relativistic stars from a specified seed isotropic star. In Section|21 
we describe the basic features of the model for a relativistic anisotropic star, 
and present the differential equations, as a first order system, governing the 
gravitational field. The algorithm that produces a new anisotropic solution 
from a seed isotropic solution is fully derived in Section El and we present 
the final result in the form of a theorem. The isotropic isothermal model 
is used to generate a new anisotropic solution in Section HI to illustrate the 
theorem. The conventional isothermal solution with a linear equation of state 
is contained in the new class of models. As a second example the interior 
Schwarzschild model is used to generate a new anisotropic solution in Section 
El The conventional interior Schwarzschild solution with constant density is 
contained in the new class of models. The details of the integration process 
for the second example is given in the Appendix for convenience. For both 
examples we study in brief the profile of the anisotropy factor. In Section (HI 
we briefly indicate that the solutions found can be related to stellar models 
used in astronomy. 
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2 Field equations 



We utilise a representation of the field equations in which only first order 
derivatives appear. This helps to simplify the integration process as pointed 
out by Chaisi and Maharaj [11] whose notation and conventions we follow. 
We take the line element for static spherically symmetric spacetimes to be 

ds^ = -e'^dt^ + eMr^ + (d^2 ^ sin^ (1) 

where i>{r) and A(r) are arbitrary functions. The energy- momentum tensor 
for isotropic matter which is not radiating has the form 

T'^'' = {fi + p)u''u^ + pg^^ (2) 

where the energy density fi, and the isotropic pressure p are measured relative 
to the four-velocity = e^^^'^S^ which is comoving. We define the mass 
function as 

1 r 

f^ij) = o / x^lJi'{x)dx (3) 



2^0 

so that M = m{R) is the total mass of a sphere of radius R. The Einstein 
field equations, with the help of (HI)-© are equivalent to the system 

e"^ = 1 4a 

r 

r(r — 2m) z/' = pr^ + 2m (4b) 

{fi + p)u' + 2p' = (4c) 

for isotropic matter distributions. We are using units where the speed of 
light (c) and the coupling constant (SttG/c^) are unity. 

The energy-momentum tensor for anisotropic matter which is not radiat- 
ing has the form 

T'^^ = {H + p)u''u^ + + tt"'' (5) 



The quantity vr"^ = \^S{r) (c"c^ — is the anisotropic stress tensor; 

the spacehke vector = is orthogonal to the fluid four-velocity 

u"" = e~'^/^(5o and |5'(r)| is the magnitude of the stress tensor. The Einstein 
field equations, with (H)) and (0), can be written as 

e"^ = 1 6a 

r 

r(r — 2m)u' = prv'^ + 2m (6b) 
4 

(ij, + Pr) u + 2p^ = — iPr-p±) (6c) 

r 

for anisotropic matter distributions. The quantities Pr and p± are the radial 
and tangential pressures respectively. The radial pressure Pr = p + 23/ -\/3 
and the tangential pressure p± = p — S/ \^ are not equal for anisotropic 
matter. Note that for isotropic matter S = and Pr = P± = P and we regain 
(jH). The magnitude S provides a measure of anisotropy. 



3 The Algorithm 

In this section we establish a procedure for generating a new anisotropic solu- 
tion of the Einstein field equations from a known isotropic solution. Consider 
the Einstein field equations (j3)) with isotropic matter distribution. Suppose 
an explicit solution to (0)) is known where 

(z/, A,m,p) = (z/Q, Ao,mo,po) (7) 
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and the functions uq, Aq, rriQ and Pq are explicitly given. Then the equations 
in (jH) are satisfied and we can write 



e 



r 

,3 



r(r — 2mo)z/Q = po'" + 2mo (8b) 
(^+^0)^^0 + 2^0 = (8c) 

The equations in (jH)) correspond to an isotropic relativistic sphere. Now 
consider the Einstein field equations © with anisotropic matter distribution. 
We seek an explicit solution to the system (jHl). To this end we propose the 
possible solution 

(z/, A, m, pr, pi_) = (^uo + p{r), Aq, mo, po + a{r), Po - ^a{r)^ (9) 

where (z/q, -^0,^0,^0) are given by ((Tj), a and /? are arbitrary functions, and 
we have set a = 2S/\/3 for convenience. The equations Q apply to the 
anisotropic relativistic sphere. 

With the specified form ©, the system © becomes 

e-^o = I 

r 

r{r — 2mo)i'Q + r{r — 2mo)P' = pqt"^ + ar^ + 2mQ (10b) 

/2m' \ 6 

i^+Po + aUiy'Q + P') + 2p'o + 2a' = --a (10c) 

The systems (jH} and (|Tn|l imply that 

(r-2mo)/5' = ar^ (11a) 
i-^+PoU' + a{u'^ + (3') + 2a' = -\ (lib) 

which is a coupled system of first order equations. It remains to integrate 
(jllaj] and obtain a and (3. 
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We can write (lllbl) as 



[^+Po] ^ + u', + [3' + 2- = — 

\ J r — zniQ a r 



with the help of ()lla|) . This differential equation can be integrated to give 
a = ^exp|-^(J„ + z/o + /3)| (12) 

where is a constant, and we have set 

2m'r. \ , , , 

-+Po] T^dr (13) 



J r — 2mo 

From (jllaj) and (fT^ we generate the nonlinear differential equation in j3: 

This nonlinear first order differential equation is integrable and we generate 
the result 

(3 = 21n|^/^ + 4 (14) 



^2 

where £ is a constant of integration, and we have set 

exp {-| {la + Z/q)} 



dr (15) 

r (r — 2mo) 

Thus the anisotropic field equations (jH)) have been integrated to generate 
a new exact solution. Note that with = 0(= a) this algorithm regains 
the isotropic solution ((7j). When 7^ then the solution is necessarily 
anisotropic. The integrations in 1^ and Ip in (fT^ and (fT^ respectively can 
be performed explicitly as z/q, po and mo are specified in the isotropic solution 
functions in (jZj). We make two observations relating to the constants of 
integration. These could have been introduced at a later stage as the solution 
is expressed in terms of indefinite integrals. If = then the parameter i has 
to be greater than zero and can be removed by rescaling the time coordinate. 



We summarise our result in terms of the following theorem: 



Theorem A: If (z^o, ^o, ^o, Po) is a given isotropic solution of the Einstein 
field equations then (z/, X,m,pr,p±)= {i^o + /5(r), Aq, mo,po + a{r),po — |a;(r)) 
is a new anisotropic solution where 



k 

a = — exp 



2m'n 



r — 2mo 



dr + uo + p 



2 In 



2m' 



Po 



r—2mo 



dr + uo 



)} 



r (r — 2mo) 



dr + 



and k and £ are constants. 



4 Anisotropic isothermal spheres 



As a first example we illustrate the applicability of Theorem A by generating 
anisotropic isothermal spheres. The line element for the isothermal model 
B is 

ds2 = -r iT^dt^ ^ ^ ^ ^^^^ ^ ^^^^ ^^^^ 

where c is a constant. The relevant isotropic functions for ()lfj|l are 



4c f 4c 

(z/Q, Ao,mo,Po) = ( — — Inr, ln<^ 1 + 



1 + c"""' ■ (l + c)2 

2cr 1 4c^ 



4c+(l + c)2' r2 4c+(l + c)2 



(17) 



The energy density function, that generates the mass function uiq, has 
the form 

4c 1 

^0 ~ ~A A2~2 

4c + (1 + c)^ 



From (fTHj) and po in (fTTj) we observe that 



Po = c/io (19) 

which is a hnear barotropic equation of state. Isothermal spheres with /i oc 
and the equation of state ()19|1 arise in both Newtonian and relativistic 



stars 
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12| . They have a long history in astrophysics as an equilibrium 
approximation to more complicated systems which are close to a dynamically 
relaxed state |l^. 

With the isotropic functions (fTTj) it is possible to evaluate the integrals 
la and If3 in (fT^ and (fT^. Then we generate the expressions 

, k 4C+(1 + Cf _l+5c\~^ 3+7c 

a = k\i h -^r r 

V 2(l + c)(l + 5c) J 

n oi ^ 4c+(l + c)2 _i+5c 

3 = 2\n<i —r 1+'= 

^ \ 2(l + c)(l + 5c) 

Hence the new line element has the form 



4c 



k4:C+(l + cV _l + 5£\\ 9 /_ 4c 



V 2(l + c)(l + 5c) ; V (l + c)V 

(d02 ^ gjj^2 ^^^2^ (20) 

and the matter variables have the explicit form 



Pr 



2cr 

4c + (1 + c) 
1 4c^ 
r2 4c + (1 + c) 



-1 

3+7c 



, , „ A; 4c+ (1 + c)^ 1+5C , „, ^ ^ ^ 

+k{i--- ^- ^r~~ r"^ 21b 

^ 2(l + c)(l + 5c) ' ^ ' 



Ac 



,2 



r2 4c+ (1 + c)2 

-r i+t: r 1+'= (21c) 



k f „ k 4C+(1 + C)^ l+5c \ 3+7c 



2 V 2 (1 + c)(l + 5c) 
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The isotropic isothermal sphere model (fT^ generates the anisotropic isother- 
mal sphere model ()20|) . With the parameter values k = 0, i = 1, we regain 
the conventional isothermal sphere. 



The degree of anisotropy is 

2^ V 2(l + c)(l + 5c) ' ^ ' 



A graph of the anisotropy factor (j22j) was plotted with the help of Mathemat- 
ica jl^. This is given in Figure H for the particular values of the parameters 
shown. The anisotropy factor S is plotted against the radial distance on the 
interval < r < 1. It is worth noting from this graph that the anisotropy has 
the feature that it is a monotonically decreasing function as r approaches the 
boundary subject to a particular choice of parameters. There is a singularity 
at r = which S shares with the other dynamical and metric functions. How- 
ever there are other choices of parameters that could be made such that S is 
a monotonically increasing function if the physics of the problem demanded 
such behavior. An increasing profile for S would be physically relevant for 



boson stars as pointed out by Dev and Gleiser p|; however note that their 
analysis was performed for a constant energy density. The simple behavior 
of S reflected in this graph indicates that a full physical investigation of this 
solution is possible; we will perform this investigation in future work. 



5 Anisotropic Schwarzschild spheres 

As a second example we illustrate the utility of Theorem A by generat- 
ing anisotropic Schwarzschild spheres. The line element for the interior 
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Figure 1: Decreasing S{r) for anisotropic isothermal sphere; c = 0.1, k 
0.005, and £ = 10 



3 if 



Schwarzschild model LLiJl is 



„2 \ -1 



(d^2 ^ gjj^2 ^^^2^ 



21n<!A-5Wl--^ , -lMl-;^r, 



r 

ZR2' 



1 /^-35Vl-i. 



A-fi./l 



ij2 



(23) 



where A and B are constants. The relevant isotropic functions for ()23|) are 

,2 

(i/Q, Ao,mo,Po) 



(24) 



The energy density function, that generates the mass function mg, has 
the form /iq = We therefore have 



constant 



(25) 



which essentially replaces the equation of state (jTHl) of Section HI The interior 
of dense neutron stars and superdense relativistic stars are of near uniform 
density 13 ■ Consequently the assurnptio n (|25ll of uniform energy density 
is often made in the modelling process jlL loL Il9l|. 
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The integral in (fT^ is easily integrated and we obtain 



(26) 



However the integral Ip takes the form 



rdr 



(27) 



The integrand in ()27|) is complicated; however the integration can be com- 
pleted and Ij3 can be written purely in terms of elementary functions. The 
details of the integration are given in the Appendix. 

We need to investigate two cases: A ^ B and A = B. We observe from 
the Appendix that la and J/j can be written in closed form for these two 
cases. Two line elements arise and we represent our solutions as follows: 

Case I: A^ B 



In this case the line element is 

^ 2 



A-B\ I 



2R{A^-B^y \ VB^ - A^ 



tanh-iJ ^^^^^'^''^^'^''''^^ 



- A^ 



^{2AB+{A^ + B^) 



+ 



(28) 
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and the matter variables become 

^3 



m 



Pr 



(29a) 



1 p-3i?v/l- 7^ 



1 



R2 



A-B^|l-^ 



rHA-B^|l-^ 



k 



6Afi2 ) (A + 5)tan{isin-^^} 



r 



(A^ + E^) Wl- 



i?2 



+ 



-1 



(29b) 



1 ^^-35^1-^ 



k^/l-§ 



/J2 



QAB'^ 



2r3 ( A-B^l 

_^j{A + B) tan{isin-i^} 



2R{A^- BY \ - 



tanh 



- A2 



-^|2A5+(A^ + 5^)Vl-^ 



-1' 



-1 -1 



With the parameter values k = 0, 
Schwarzschild sphere 



(29c) 



1 we regain the original interior 
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The degree of anisotropy is 



S 



/?2 



2r3 {A-B^/l-j^ 



X 



k 



6AB^ 



tanh ^ 



(A + B) tan{isin-i j^} 
- A^ 



^{2AB+{A^ + B^) 



-1 -1 



B^r 



A-B\ I 



i?2 



+ 



(30) 



The graph of the anisotropy factor was plotted with the assistance of 
Mathematica Q| . This is shown in Figure El for the particular values of the 
parameters. The interval for the plot of S against r is < r < 1. The 
quantity is a monotonically decreasing function. Subject to the choice of 
the parameters, the anisotropy can be constructed such that it is a mono- 
tonically decreasing function as r approaches the boundary. The function 
S vanishes at the boundary. Other choices of the parameters A, B, k and 
i may generate different behavior for S. The singularity at the center does 
not seem to be 'removable' by any choice of parameters. A comparison of 
this case with Dev and Gleiser |lj] is not possible as the anisotropy factor 
vanishes at the boundary r = R. This will be investigated further and a 
detailed analysis of the anisotropy factor S and the dynamical variables for 
this anisotropic solution will be pursued in the future. 
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0.2 0.4 0.6 



Figure 2: Decreasing S{r) for anisotropic Schwarzschild sphere (^4 7^ B); 
A^1,B^ -10, A; = 2, ^ = 1 and i? = 1. 

Case 11: A^B 

For this case the hne element is 



X 



X 



80R 



i?2 



1 . _i r I \l . _^ r 
CSC" < - sm — > sec < - sm — 



R 



R 



'-w)y'-R^-'['-w) 



\ -1 
r \ 



+ ( 1 - — J dr^ + {de^ + sin^ ^d0^) 



(31) 
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and the matter variables are given by 
m 



(32a) 



Pr 




1 . _i r I J 1 . _i r 



2 1 



2r2 



(32b) 



CSC < - sm — > sec < - sm — 

80R 12 R\ \2 R 



X 



2r2 



2r2 



(32c) 



With the parameter values A; = 0, £ = 1, we regain the original interior 
Schwarzschild sphere (I23|l with A = B. 



The degree of anisotropy is 



S 



2r3 1 



1 



CSC < - sm — > sec < - sm — 

80R 12 R\ \2 R 



2r2\ 



(33) 



The graph of the anisotropy factor (j33|l was plotted with the assistance of 
Mathematica [15] . This is shown in Figure El for the particular values of 
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the parameters. The interval for the plot of 5* against r is < r < 1. 
The quantity S* is a monotonically decreasing function. The behavior of 5* 
is similar to the case A ^ B given in Figure El However observe that the 
behavior in Figure El is more restricted for Case 11 as A and B are fixed. 

2 

1 . 5 

1 

CO 

0.5 



0.2 0.4 0.6 0.8 1 

r 

Figure 3: Decreasing S{r) for anisotropic Schwarzschild sphere {A = B); 
k = 0.00006, i = 100 and i? = 1 



6 Discussion 



It is remarkable that the integrations in Theorem A can be completed in 
closed form in spite of the complicated nature of the integrands in a and /3, 
and we obtain expressions in terms of elementary functions. It is possible to 
demonstrate that the solutions found can be used to discuss the structure 
of neutron stars and quasi-stellar objects. Now consider a neutron star of 
radius 10 kilometres and surface density of 2x 10^^ gcm~^ in c.g.s units. Then 
the parameters relating to the surface redshift z = {1 — 2M/R) — 1, and 
the mass M in terms of the solar masses Mpt can be calculated. We choose 
values so tKat co.panson w.th Gokh.oo a.d Meh.a R is facil.tated, Fo. 
the choice c = 0.1, we obtain z = 0.154 and M = 0.520Mq. The choice 
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c = 0.9, generates the values z = 0.413 and M = 1.40Mq. Clearly other 
values for z and M can be generated by making different choices for c and 
the surface density. The range of values obtainable is consistent with the 
results of Gokhroo and Mehra Ql- Hence our solutions yield values for 
surface redshifts and masses that correspond to realistic stellar sources such 
as Her X-1 and Vela X-1. 
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Appendix 



The integral J/3 has the form 



rdr 



Note that two cases arise in the integration process: A B and A = B. 
Case I: A^ B 

To carry out the integration in (jTfj) for A ^ B we make the substitution 

r 

sint? = — 
R 
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so that //3 becomes 



1 



RJ sin2 7?(A-5cost9)2 
1 



- {2AB + A^ cos d + B"^ cos d) esc d 



5^ sin -i? 



A - Bcosd 

in terms of "(9. In terms of the original radial coordinate r we have 
1 



'/3 



6A52 , 1 f (A + 5)tan{isin-i^} 

tanh"^ '-2 ^ J 



Vfi2 _ ^2 



— (2AB + (A^ + B^) cos |sin-i 



R{A- Bcos{sm~^ j^]) 



1 



^tanh-^J(^ + ^)t^^{i^i^"'i} 



- A2 



2 



Then ((121) and (HH) give 



i?2 



rHA-B^|l-^ 



X 



2R{A^-B^y \ V52 - A2 



tanh-^J (^ + ^)tan{|sin-^^} 



V52 - A2 



2 



+ 
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R \ V 
for the functions a and /5 



Case II: A = B 



To carry out the integration in (P7j) with A = we again make the 
substitution 

r 



sin ^9 

R 



to have 



rdr 



-R 7 sin^ (1 — cos'i?) 

1 {} 

csc^ — sec — (5 cos iD — A cos 2'i9 + cos S-i?) 



80i? 2 2 

The terms cos 2d and cos "id can be simphfied with basic trigonometric iden- 
tities and J/3, in terms of the original radial coordinate r, becomes 



1 5 / 1 • -1 ^ 1 f 1 . -1 ^ 



,2 
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Then (HH) and (HH) give 




for the functions a and /3. 
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